Proton Fraction in Neutron Stars 
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The proton fraction in /3-stable neutron stars is investigated within the framework 
of the Skyrme-Hartree-Fock theory using the extended Skyrme effective interaction 
for the first time. The calculated results show that the proton fraction disappears at 
high density, which implies that the pure neutron matter may exist in the interior of 
neutron stars. The incompressibility of the nuclear equation of state is shown to be 
more important to determine the proton fraction. Meanwhile, it is indicated that the 
addition of muons in neutron stars will change the proton fraction. It is also found 
that the higher-order terms of the nuclear symmetry energy have obvious effects on the 
proton fraction and the parabolic law of the nuclear symmetry energy is not enough 
to determine the proton fraction. 

PACS: 97.60.Jd, 26.60,+c 

The neutron star serves as a cosmic laboratory to study the nuclear equation of state (EOS) at high 
density. It has been known that the mass, radius, and density profile of neutron star are determined by 
the EOS of asymmetric nuclear matter, and the chemical composition, particularly the proton fraction 
is determined by the nuclear symmetry energy 1 . Many theoretical methods, such as the Brueckner 
approach 2 , variational many body (VMB) theory 3 , the relativistic mean field (RMF) theory 4 , and 
non-relativistic many body theory 5 , have been used to study the proton fraction in neutron stars 
since the proton fraction is an important physical quantity to determine the evolution of neutron 
stars, such as the cooling rate, neutrino flux, and so on 6-9 . However, little is known about the 
properties of the nuclear symmetry energy at high density 10 with the proton fraction in neutron stars 
remaining undetermined. 

In our previous work 11 ' 12 , a density, temperature and isospin dependent nuclear EOS has been 
deduced within the framework of the Skyrme-Hartree-Fock theory using the so-called extended Skyrme 
effective interaction 13 . For asymmetric nuclear matter with density p and neutron excess 5 = (p n — 
Pp)/p i the energy per nucleon at temperature T can be written as 
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where t q and — t q represent neutron and proton, respectively. Every equation is split into two equa- 
tions, which are for the neutron and proton, respectively. The to — t$ and xq — X5 arc Skyrmc potential 
parameters, A and Ci(p Ta ) represent the thermal wave length and Fermi-Dirac integral. The Taylor ex- 
pansion of Eq. (1) with respect to 5 at temperature T = (hereafter, we define e(p, S) = e(p, 6, T = 0)) 
is as follows 



e( P , S) = e(p, S = 0) + C{ P )S 2 + D(p)S 4 + 0(6 6 ), 



(2) 



where C is the so-called coefficient of symmetry energy strength. 

The chemical composition of the neutron star is determined by the requirement of charge neutrality 
and equilibrium with respect to the weak interaction ( /3-stable matter ). From Eq. (1) or (2) we can 
calculate the proton fraction, x p = (1 — S)/2, for /3-stable matter as found in the interior of neutron 
stars. Equilibrium for the reactions n — > p + e~ + v e and p + e~ — > n + v e requires that 

Hn = Pp + Pe, (3) 

where pi = de/dNi (i — n, p, e and p) is the chemical potential. Consequently, one obtains 
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For relativistic degenerate electrons, 
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with h = c = 1 and x p = x e because of charge neutrality. From Eqs. (1), (3)-(5), it is easy to find 
numerically the proton fraction as a function of density p. 

Just above nuclear matter density at which p e exceeds the muon mass m M , the reactions e~ — ► 

> p + p~ may be energetically allowed such that both 
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electrons and muons are present in /3-stable matter. This alters the /3-stability condition to 
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with x p = x e + x^ 

The bulk properties of symmetry nuclear matter at saturation density, such as the density p , the 
binding energy per nuclcon e 0j the incompressibility K, the effective mass m*/m, and the symmetry 
energy are summarized in Table 1 for the three parameter sets SII, SHI, and SkI5. From Table 1, one 
can see that the difference for incompressibility between SII and SkI5 is very large, about 85 MeV 
at saturation density, while they almost have the same m*/m=0.58. Meanwhile, the difference for 
m*/m is very large between SIII and SII, about 0.18 at saturation density, while the difference for 
incompressibility between them is relatively small, about 14 MeV. In calculations, we note that the 
binding energy per nucleon e and the nuclear symmetry energy C exhibit completely different behaviors 
with the increase of the baryon density for different Skyrme parameter sets. These differences come 
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from the difference of the nuclear EOSs. Therefore, one may investigate the proton fraction for 
different nuclear EOSs. 

Under the parabolic approximation (up to the -S 2 term in Eq. (2)), the chemical compositions of 
neutron stars in both the electrons only and electrons plus muons cases as a function of the baryon 
density are shown in Figs. 1 (a)-l (c) using the Skyrme parameter sets SII, SIII, and SkI5, respectively. 
In order to ensure validity of the proton fraction above, we show the squared ratio of sound speed 
to light speed, (s/c) 2 = dP/dE (where E is the total energy density), in Fig. 1 (d). The causality 
condition requires (s/c) 2 < 1. One notes that for SII and SkI5, the squared ratio (s/c) 2 will make 
the causality condition fail to be satisfied at baryon density larger than 1.35 fm~ 3 and 1.1 fm~ 3 , 
respectively, about 8 times larger than the saturation density for symmetry nuclear matter. For SIII, 
up to baryon density 1.5 fm~ 3 , the causality condition can be still satisfied. From Figs. 1 (a)-l (c), 
one can find that the addition of muons increases significantly the proton fraction over a broad range 
of density. In case of SIII, the muon fraction is very small and muons only exist at about 0.16 fm~ 3 , 
which is due to the fact that p e is too small to exceed the muon mass m M . The value of the proton 
fraction is an important ingredient to model the thermal evolution of a neutron stars. In fact, if the 
proton fraction in the interior of star is above a critical value x Urca , the so-called direct Urea processes 
n — > p + e~ + V e and p + e~ — > n + v e can occur 14 . If they occur, neutrino emission and neutron 
star cooling rate are increased by a large factor compared to the standard cooling scenario 14 . The 
critical proton fraction has been estimated to be in the range of x Urca = 11 — 15%. It is indicated 
from Fig. 1 that the proton fraction for SII and SIII cannot reach this range. However, the proton 
fraction reaches x Urca = 11% at baryon density 0.2 fm ~ 3 for SkI5. This means that one can get a 
reasonable proton fraction for a softer nuclear EOS, such as K = 255 MeV for case of SkI5. While 
for stiffer ones, such as SII (K = 341 MeV) and SIII (K = 355 MeV), one might not get reasonable 
proton fraction although they may have a very different effective mass (m* jra = 0.580 and 0.763 for 
SII and SIII, respectively). It seems that the incompressibility is more important to determine the 
proton fraction. 

One of the interesting features of Fig. 1 is the disappearance of the proton fraction at high density, 
namely, p — 0.48 and p = 0.32 fm~ 3 for SII and SIII, respectively. This phenomenon may be due 
primarily to the greater short-range repulsion in isospin singlet nucleon pairs compared to isospin 
triplet nucleon pairs 15 . At high density this short-range repulsion must dominate and pure neutron 
matter is favored. At intermediate densities the strong isospin singlet tensor potential and correlation 
serve to keep the isospin singlet nucleon pairs, and thus symmetry nuclear matter, more attractive 
than pure neutron matter. These features imply that the pure neutron matter may exist in the interior 
of neutron stars where the density is high enough. This phenomenon has also been observed in the 
calculations of the VMB theory. For parameter set SkI5, the disappearance of the proton fraction is 
not found up to 1.5 fm~ 3 . 

Many theoretical studies have shown that the empirical parabolic law of the single particle binding 
energy is well satisfied within a broad range of density and neutron excess 5, namely, in Eq. (2) the 
S 2 term is significant, but the higher-order terms are very small. However, it is still unknown how 
the higher-order terms affect the proton fraction at high density. Fig. 2 is the same as Fig. 1 except 
for including all the higher-order terms. From the comparison between Figs. 1 and 2, one can find 
that the addition of the high-order terms increases significantly the fraction of different constituents, 
p, e~ and pT , over a broad range of density for parameter sets SII and SIII. It implies that the 
higher-order terms of the nuclear symmetry energy have obvious effects on the chemical composition 
and the parabolic approximation of the nuclear symmetry energy may result in large deviation. 
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TABLE CAPTIONS 
Table I Bulk properties of symmetry nuclear matter at saturation point. 



Parameter 


SII 


SIII 


SkI5 


po (fm- J ) 


0.1484 


0.1453 


0.1558 


£o (MeV/nucleon) 


-15.99 


-15.85 


-15.85 


K (MeV) 


341.38 


355.35 


255.78 


in* /m 


0.580 


0.763 


0.579 


C (MeV) 


34.26 


28.25 


36.75 



FIGURE CAPTIONS 

Fig. 1 With parabolic approximation, the compositions of neutron stars in both the electrons only 
and electrons plus muons cases as a function of the baryon density using the Skyrme parameter 
sets SII (a), SIII (b), SkI5 (c). Meanwhile, the squared ratio of sound speed to the light speed 
(s/c) 2 as a function of the baryon density for these three parameter sets is displayed in (d) in 
electrons plus muons case. In (b), the muon fraction is very small so that the solid line, dotted 
line and dash-dotted line coincide with each other. 

Fig. 2 Same as Fig. 1 but including all higher-order terms instead of parabolic approximation for 
nuclear symmetry energy. 
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